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Abstract 

As an applicaton of sonification, a simulation of the early universe was de- 
veloped to portray a phase transition that occurred shortly after the Big Bang. 
The Standard Model of particle physics postulates that a hypothetical particle, the 
Higgs boson, is responsible for the breaking of the symmetry between the electro- 
magnetic force and the weak force. This phase transition may have been respon- 
sible for triggering Baryogenesis, the generation of an abundance of matter over 
anti-matter. This hypothesis is known as Electroweak Baryogenesis. In this simu- 
lation, aspects of bubble nucleation in Standard Model Electroweak Baryogenesis 
were examined and modeled using Mathematica, and sonified using SuperCol- 
lider3. The resulting simulation, which has been used for pedagogical purposes 
by one of the authors, suggests interesting possibilities for the integration of sci- 
ence and aesthetics as well as auditory perception. The sonification component in 
particular also had the unexpected benefit of being useful in debugging the Math- 
ematica code. 

1 Introduction and Motivation 
1.1 Motivation for Sonification 

Questions regarding theoretical particle physics and cosmology are of increasing in- 
terest to the scientific community due to the current operation of the Large Hadron 
Collider. On the experimental end, sonification has been explored as a tool for particle 
detection [1 1. On the theoretical end, sonification has a history of application through 
simulations of rigorous and complex theoretical models 12) Q. Perhaps the most im- 
portant benefit of sonification to the scientific community is an informed and interested 
public. Sonifications, visualizations, and simulations provide a point for public ac- 
cess to the otherwise challenging aspects of contemporary physics. Additionally, the 
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stochastic and probabilistic nature of this process lends itself to discussions of aesthet- 
ics in the sciences Pfl . 

This simulation was developed to portray the universe as it underwent a fundamen- 
tal change to a lower-energy state about 1CP 10 seconds after the Big Bang. This "phase 
transition" lasted 10 -13 of a second and left the universe with a permanent abundance 
of matter over anti-matter. The simulation illustrates this very brief but important point 
in the evolution of the universe. 

In this simulation, the phase transition has been portrayed as proceeding by bubble 
nucleation. Similar phase transitions are familiar from everyday life; as water vapor 
cools in a cloud, the water condenses into its liquid phase. Impurities or dust particles 
in the cloud serve as centers for the nucleation of bubbles, and the bubbles expand 
as the surrounding air molecules are tipped into the lower energy state, triggering a 
first-order phase transition from liquid to solid. 

For the purposes of this research, we assume that a similar first-order transition 
took place in the early universe. Our final product is a set of 13-second movie files 
(.mov) showing circular bubbles (representing areas of broken electroweak symmetry, 
the lower-energy state) appearing in a background of static (representing energy fluctu- 
ations in a hot universe). The timing and location of the bubble formation is governed 
by a probability function rather than hard-coded, so that each time the code is run a 
different distribution of bubbles result. The formation of each bubble is accompanied 
by a "chime" carrying audio information about the position of the bubble. As the bub- 
bles expand and fill the screen, there is an increase in volume and change in timbre to 
reflect that more and more of the universe has transitioned to the lower-energy state, 
producing the abundance of matter over antimatter. 



1.2 Conditions for Baryogenesis 

Although antimatter can be created in small quantities in particle accelerators, by cur- 
rent observations the visible universe is practically void of anti-matter on all scales. 
Most of the mass in the universe is baryonic. A baryon, of which protons and neutrons 
are the most common example, is a particle made up of 3 quarks. Anti-baryons are 
composed of 3 anti-quarks. The asymmetry between the number of baryons per unit 
volume (n b ) and the number of anti-baryons per volume ( ) is defined as 

n b - n b 
n b + n b 

In practice, baryons will annihilate with anti-baryons until the supply of anti- 
baryons is exhausted, producing photons: 

B + B ->• photons. (2) 
The baryon asymmetry is given experimentally by the baryon-to-photon ratio: 

v = ^ = 6 • 1(T 10 , (3) 

where n 7 is the number density of cosmic microwave background photons, assumed to 
have been produced by the annihilation process in eq. [2] Thus, for every billion baryons 
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originally in the universe, there must have been a deficit of about one anti-baryon. This 
matter- antimatter asymmetry remains one of the most important unanswered questions 
in cosmology. 

Historically, the abundance of matter was attributed to the initial conditions of the 
Big Bang. In 1967, Andrei Sakharov [5] proposed three conditions that would allow 
for Baryogenesis, the production of more matter than anti-matter: 

A) baryon number violation in the laws of nature, 

B) CP- Violation, and 

C) a departure from thermal equilibrium. (4) 

Baryon number is defined as the number of baryons minus anti-baryons in the uni- 
verse. Sakharov's first condition requires that there be some process that can gener- 
ate an excess of baryons. CP- Violation, or charge-parity violation, occurs when left- 
handed particles and right-handed antiparticles do not behave in the same way. If 
there is not enough CP- Violation, the anti-matter process that generates an excess of 
anti-baryons would occur at the same rate as the process that generates an excess of 
baryons; thus no imbalance could be formed. Finally, if the universe were in thermal 
equilibrium, the reverse reactions would occur at the same rate. Hence, a departure 
from thermal equilibrium is necessary to achieve the matter anti-matter imbalance. All 
three of these conditions may be applicable to the Electroweak Phase Transition, lead- 
ing to speculation that the breaking of the electroweak symmetry may have been the 
source of the matter-antimatter asymmetry. 

1.3 The Electroweak Phase Transition & Baryogenesis 

When the universe was very young, it was highly energetic. At these high energies, 
the electromagnetic and the weak nuclear force were symmetric and mediated by four 
massless gauge bosons. The Standard Model of particle physics hypothesizes the ex- 
istence of an additional boson (the Higgs particle); as the universe cooled, the Higgs 
field would have broken the electroweak symmetry as it transitioned to its vacuum 
state, provided that its vacuum energy was non-zero. (The actual vacuum expectation 
value that fits current data is 246 GeV.) This breaking of electroweak symmetry left 
the photon massless and gave mass to the W ± and Z bosons of the weak force. This 
moment is of particular interest to Baryogenesis. The Higgs field does not acquire its 
vacuum energy uniformly throughout space; in the scenario of electroweak baryogene- 
sis, regions of space where the Higgs field transitions to its vacuum state form bubbles 
of broken symmetry, which would expand or shrink depending on their radii. This pro- 
cess is called bubble nucleation and it is within the bubble "walls" that Baryogenesis 
could occur. 

All three of Sakharov's condtions may be present in the Electroweak Phase Tran- 
sition. Baryon number violation occurs in standard electroweak theory; although it 
is suppressed at low temperatures, baryon number violation becomes stronger at high 
temperatures [6 1. The following toy example illustrates how Electroweak Baryogenesis 
could proceed: Suppose the process to generate an excess of baryons is X — > Y + B. If 
there is enough CP- Violation in the model then this process may occur more often than 
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its anti-matter equivalent process X — > Y + B. If the phase transition is significantly 
first order, or extremely abrupt, the region of rapidly changing Higgs field would be 
out of thermal equilibrium. Hence the reverse processes Y + B — > X and Y + B — » X 
would not occur at the same rate as the forward processes, perpetuating the imbalance. 
The excess of baryons would be swept into the bubble and acquire mass. Since the new 
conditions are drastically different inside the bubble, the excess matter would persist 
and baryon number would be violated. Theoretical studies of baryon number violating 
processes have shown that such a process is possible within the walls 0. 

In this simulation, we assume, for simplicity, that the breaking of electroweak sym- 
metry is governed by the Higgs content of the Standard Model (a single Higgs parti- 
cle). This is not the most realistic electroweak baryogenesis scenario; for one thing, 
the Standard Model does not provide enough theoretical CP- Violation for Electroweak 
Baryogenesis to generate the observed excess of matter. Furthermore, the mass of the 
Higgs boson would have to be less than roughly 40 GeV to create a sufficiently large 
departure from thermal equilibrium (5). This mass is improbable because particle ac- 
celerators would have already discovered a Higgs particle this light; an electroweak 
phase transition in line with current experimental bounds on the Higgs mass would be 
at best only weakly first order. Other models such as Supersymmetry or a Two-Higgs 
Doublet Model may allow for more CP- Violation and possibly a higher Higgs mass. 
Thus, Electroweak Baryogenesis is not ruled out in those models. 

This paper provides a model of bubble nucleation in Electroweak Baryogenesis 
assuming a Standard Model Higgs sector and a Higgs mass of 35 GeV. Although this 
mass is unrealistic according to current experimental bounds, it is a useful starting 
point for modeling Electroweak Baryogenesis. The mathematical formalism will be 
simplified in favor of a more rigorous explanation of the sonification techniques. For 
the full mathematical formalism see ref. (9). 

2 Effective Potential, Critical Radius, & Wall Velocity 

The evolution of the Higgs field with time (or equivalently, temperature) is governed by 
its effective potential (see Fig. [TJ. The equations used to define the effective potential 
in this work were taken from |8 1. To create the simulation, one must derive values for 
the critical radius and wall velocity of the bubbles. For a Higgs mass of 35 GeV, the 
critical temperature was found to be T c = 71 .4 GeV. This is the temperature at which 
a second minimum forms in the potential. Below this temperature, the state of broken 
electroweak symmetry (a non-zero vacuum expectation value for the Higgs field) is 
energetically favorable. 

Based on this effective potential, the critical radius was derived. The critical radius 
(which depends on temperature) is the minimum size of bubble that will allow the 
bubble to survive and grow; bubbles that are nucleated with subcritical radii do not gain 
enough energy to overcome surface tension, and thus shrink and disappear. Equations 
used for critical radius in this simulation were taken from ref. ifTTl . 

The derivation for the wall velocity was found in lfT2l . Discrete values of wall 
velocities for a corresponding Higgs mass from [ 12 1 were interpolated in Mathematica 
in order to find a value for a Higgs mass of 35 GeV. The wall velocity was approximated 
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Figure 1: The effective potential for a Higgs mass of 35 GeV at 71.4 GeV tempera- 
ture. As the universe cools to the "critical temperature" T c , the Higgs potential forms a 
second minimum. The deepening of this minimum (shown here at 60GeV) as the cool- 
ing continues forms a "true vacuum" (global minimum). Energy fluctuations allow the 
Higgs field to transition into the true vacuum, forming a bubble of broken electroweak 
symmetry. 



to be V w = 0.375 c, where c is the speed of light. 



3 Probability Functions 

Now that the critical radius and wall velocity have been described, the probability of 
bubble formation at a particular time, and the initial radius of a bubble are also needed 
for the simulation. 

3.1 Determining Initial Radius 

To determine the initial size of a bubble, a generic approximation of bubble nucle- 
ation was applied to the Electroweak Phase Transition. The probability function, P r , is 
defined in |[T3l and applied to the simulation. 

When the function is plotted, the area under the curve represents the probability that 
a bubble will be formed within the corresponding domain. We note that the function 
is also temperature dependent so the distribution of bubbles will not be uniform with 
time. A higher probability is given for a bubble to shrink near the start of the phase 
transition. However, as the temperature cools the function gives equal probabilities for 
radii sizes. 
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3.2 Determining Formation Time 



The electroweak theory determines the probability that a bubble will form per unit time 
and per unit volume. As the universe evolves in time, it also cools, so time is related to 
temperature. An approximation for the probability of formation with temperature, P t , 
was found in lfl4ll and applied to the simulation. 

This probability function works the same way as P r in the way of determining the 
probability of a particular domain. By examining Fig. [2] it can be seen that bubbles are 
more likely to form near the end of the phase transition. 

Probability 




Coefficent of Tc 



0.9885 0.9890 0.9895 0.9900 

Figure 2: The probability of bubble formation vs. temperature. Values before 0.9882 
T c were complex. 



4 Programing the Simulation 

The goal of this simulation was to demonstrate how these bubbles might form and grow 
over a certain area of space and time. 

4.1 Scaling Space and Time 

In order to have a physical simulation the correct scaling for space and time had to be 
determined. For scaling time and temperature, the domain of the real values for P t was 
used. These ranged from 0.99 T c - 0.9882 T c , This was then converted to a 10 second 
time scale. 

A way to establish the spacing of the bubbles was also needed for the simulation. 
From [14 1, the equation r s will give the average spacing between a bubble and the 
nearest bubble. When thirty bubbles were programmed to form within a 300 x 300 unit 
in 2 dimensional space, the average spacing was found to converge to 29.7 units. 
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4.2 Generating Parameters 

The next step was to program how to determine the parameters for a given bubble. 
The X and Y coordinates for the center of each bubble were randomly chosen with no 
weights. This proved to be the simplest method and since thirty bubbles were used, an 
even distribution often followed. The formation time of a given bubble was weighted 
by Pt and the initial radius by P r . Next, an algorithm to reject bubbles that would form 
overlapping with other bubbles was written. 

After all of the parameters were calculated, they were exported to be simulated. 
The parameter set next was manipulated to the necessary format for sonification and 
then exported in the form of .csv files for the sonification process. 

5 Sonification 

A selection of 15 simulations have been made available on the video-sharing site 
"YouTube," and can be found under the keyword 'bubble phase transition.' Screen- 
shots from the 9th simulation are displayed in Fig. [3] 

t w 4(s) t w 7(a) t « 8(a) t w 9(a) 

Figure 3: Four screenshots from the 9th simulation. Before t ps 4 seconds, there has 
been no bubble nucleation. After t « 10 seconds, the phase transition is complete. 

In order to enhance the simulation, a sonification was created using the program 
SuperCollider 3JT3J. Two synthesis definitions were employed, one to represent the 
formation of the bubbles over time, the other to represent the amount of space that still 
had Electroweak Symmetry. The sonification techniques employed were judged to be 
effective based upon brief discussions of auditory perception with Dr. Neuhoff and Dr. 
Walker in 2009, various papers from the ICAD proceedings, and the current level of 
experience with SuperCollider. However, insights from the international community 
would be beneficial to ongoing research goals. 

5.1 Sonification of Bubble Formation 

The sonification of bubble formation was modeled upon three known quantities, 

A) the time that the bubble first appeared, 

B) its horizontal (X) position, and 

C) its vertical (Y) position, (5) 

as generated using the previously described Mathematica program. Mathematica wrote 
these values into a .csv file which was imported into SuperCollider. Using A from 
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two additional values were derived, 

D) time until the end of the simulation, and 

E) time between bubble events. (6) 

An example data set imported into SuperCollider is displayed in Table [T] 

Table 1: An example data set imported into SuperCollider for a 13 second simulation. 
A, B, C, D and E refer to the titles in Q and |6]). The values for A, D, E are in seconds 
and C, D are the relative X and Y positions on a 300 x 300 visual grid. 



Bubble # 


A 


B 


C 


D 


E 


1 


5.6 


23.44 


73.45 


7.4 


0.53 


2 


6.13 


254.97 


91.67 


6.87 


1.79 


3 


7.92 


193.55 


40.59 


5.08 


0.08 


4 


8 


172.75 


151.04 


5 


0.06 


5 


8.06 


89.96 


183.79 


4.94 


0.71 


6 


8.77 


161.62 


65.08 


4.23 






As discussed in previous sections, each bubble that formed had an initial radius r 
that would shrink or grow depending upon its magnitude. Values for r that were below 
some critical radius r c would briefly appear and then disappear, whereas bubbles with 
r > r c lasted until the end of the simulation. The bubbles which quickly shrank were 
not sonified because of lack of time and because a sufficient technique for sonifica- 
tion was not determined. However, because instances of shrinking bubbles (r < r c ) 
were perceptually insignificant compared to the bubbles which expanded (r > r c ), our 
choice was not thought to detract from the overall simulation. 

5.1.1 Spatialization, Pitch, Timbre 

To minimize perceptual interference between parameters and maximize aesthetic ap- 
peal, Spatialization, Timbre, Intensity and Pitch were the auditory parameters chosen 
for representation. 

For every bubble with r > r c there was one corresponding sound event such that 
the bubble's X-position corresponded to horizontal spatialization, its Y-position corre- 
sponded to pitch (high Y-position — » high pitch), and the increasing radius corresponds 
to a timbre and intensity which changes in a way intended to simulate the feeling of 
sudden appearance followed by expansion. 

In SuperCollider, spatialization was controlled using 'Pan2', a two channel equal 
power panner function. Pitch was controlled using 'Midicps', a command which con- 
verts MIDI notes to their corresponding frequency. The MIDI range used in the present 
sonification was 40-60 or alternatively, E2 to C4 in musical notation. Timbre was con- 
trolled using 'Formant,' a command which generates a set of harmonics around a for- 
mant frequency given a certain fundamental frequency. Both the formant frequency 
and the pulse width frequency increased from 100-1000 Hertz along an exponential 
curve ('XLine') in the time frame of D-l seconds. D-l refers to the first element of the 
D column in Table[T] Only the first bubble reaches the 1000 Hz level. 



8 



So that the sonification would outlast the Electroweak Phase Transition, an extra 
three seconds were added to the simulation. This additional time created a fuller audi- 
tory experience. Because the data was not generated for a 13 second visualization, only 
visualizations in which the entire visual grid had undergone transition after 10 seconds 
were deemed physical. The majority of the visualizations did not complete the phase 
transition within 10 seconds, consistent with fl4l . 

5.1.2 Intensity 

The intensity of a bubble forming was also varied over time such that the increasing 
bubble size would have an analogous auditory experience. In SuperCollider, the En- 
velope Generator ('EnvGen') was used to control the intensity over time. An example 
envelope is displayed in Fig. [4] The envelope was designed such that every bubble 
could be heard when it first entered, but by the end of the simulation, would not be 
as loud as the bubble which began before it. The intensity was made to increase or 
decrease cubically to reflect the fact that the volume of a bubble is proportional to r 3 . 




Figure 4: An example envelope for controlling the intensity over time of each bubble. 
This envelope corresponds to Bubble #4 of Table[T] The initial intensity of 1 decreases 
to 0.3 cubically in the 1 time frame of second. Having reached its minimum value it 
approaches its final intensity of 5/7.4 cubically in the time frame of 5-1 seconds. Note 
that the entire envelope occurs over the time frame of D-4 of Table [T] 

While designing the sonification, it was determined that bubbles with lower pitch 
(lower Y-coordinate), sounded softer than higher pitch bubbles (higher Y-coordinate). 
Though the exact relationship for perceived equal intensity across octaves was not 
known, an additional intensity function was created to compensate for this inequality. 
The function was constructed so that the amplitude of the lowest possible pitch (E2) 
was twice that of the highest possible pitch (C4), with a linear relationship between the 
endpoints. 
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5.2 Sonification of Electroweak Symmetry 



If the sonification involved only bubble formation, about the first five seconds of the 
simulation would be completely silent and featureless, somewhat at odds with the real- 
ity of a hot universe streaming with particles moving at the speed of light. To represent 
the fluctuations in energy present in this phase of the universe, a crackling noise was 
chosen for the area which had not yet undergone transition. In SuperCollider, this 
crackling noise was generated using 'Dust2', a function which generates a certain den- 
sity of random impulses from -1 to 1. The average number (density) of these random 
impulses per second was chosen to be a constant 1000. 

As bubbles are formed, the area which still has has electroweak symmetry slowly 
decreases. For this period of time, the intensity of the dust decreases proportional to the 
amount of space that still has electroweak symmetry. As in section [5T| the Envelope 
Genererator 'EnvGen' was used for its control of intensity. An example envelope can 
be viewed in Fig. [5] 




Figure 5: A screenshot from SuperCollider displaying the envelope of 'Dust2' over a 
period of 13 seconds. Beginning with an initial intensity of 1, as electroweak symmetry 
is broken the envelope decreases to a final value of when electroweak symmetry is 
completely broken. 



The envelope was created using the arrays generated in section 4.2 that determined 
whether bubbles would form inside other bubbles. Though more rigorous calculations 
could be employed to calculate this area over time, they were not attempted because 
of lack of time. Furthermore, the perceivable difference in sound between the two 
methods was not thought to be substantial enough to merit further inquiry. 



5.3 Analysis of Sonification 

The goal of the sonification was to augment a scientific visualization for a fuller and 
more engaging perceptual experience. The most important data properties were the 
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bubble's horizontal and vertical position, the amount of unbroken electroweak symme- 
try remaining, and the rate of expansion for each bubble. The horizontal and vertical 
position of the bubble centers, indicated by the spatial positioning and pitch respec- 
tively, communicates the areas of space that have undergone baryogenesis in each sim- 
ulation. This data property is complemented by the volume envelope, providing an ad- 
ditional cue as to the total amount of remaining electroweak symmetry. The perceptual 
feeling of expansion (increasing volume) communicates the dynamic characteristics of 
the bubble wall. 

The sonification created to model electroweak symmetry was thought to be effec- 
tive because it communicated the important data properties in an intuitive way. Its 
effectiveness was confirmed when errors in the Mathematica code were discovered 
through instances in which the combined visualization and sonification were not of the 
same data set. In one instance, a bubble was heard forming inside of another bubble 
(an unphysical event), which was not evident from the visual display. This discrepancy 
indicated an error in the Mathematica code that may otherwise have gone undetected. 

Several issues involving the intensity of the sound files arose. It was discovered that 
lower pitches were softer than higher pitches despite being equal in all other respects. 
The SuperCollider UGen 'AmpComp' provides basic psychoacoustic amplitude com- 
pensation and would have helped to provide more uniformity. Further progress towards 
loudness scaling for digital synthesis has been realized in |[T6ll . 

Though the envelope of Fig. |4] gave special emphasis to bubbles that had just 
formed, at some points it proved insufficient in this respect. There were instances 
in which new bubbles were not heard as clearly as others because of the intensity of 
the combined sounds of other bubbles which had previously formed. Though perhaps 
these new events are less relevant when there are many bubbles expanding simultane- 
ously, an envelope in Fig. [4] that started from with a shorter attack time would have 
increased the perception of new bubbles. 

Physically speaking, at the phase transition, the universe does not cool drastically. 
While the sonification effectively demonstrated the phase transition, the intensity of 
the dust sound is not as great as the intensity of the sounds after the transition. This 
discrepancy might be misleading and could be altered to keep a constant intensity over 
time. However, this adjustment might have lead to a decrease in perception of bubble 
formation and ultimately was not used. 

Another option for displaying the amount of electroweak symmetry over time 
would be altering the density of 'Dust2' over time instead of its intensity. Although 
this possibility was considered, it would have required the ability to do curve fitting in 
SuperCollider, a feature it does not support at this time. 

6 Conclusions 

The sonification presented was effective as a complement to visualization and as a 
tool for experiencing the dynamics of a phase transition that occurred shortly after the 
Big Bang. The overlap of sonification and visualization offered a means of finding 
and correcting bugs that could not have been discovered through visualization alone. 
The sonification created would benefit from insights and corrections from the science 
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of auditory perception. Simulations such as these have potential usefulness in public 
outreach and physics education. 

7 Future Work 

The actual process of baryogenesis occured at a rate 10 trillion times faster than that 
displayed in our simulation. If this simulation were slower or faster by a power of 10, 
humans likely find it less aesthetically pleasing and more difficult to interpret. How- 
ever, much finer adjustments can be made to this simulation. Is there an ideal rate that 
sonic events in this simulation (or any simulation) would occur and evolve? 

A broader area of research in the intersection between science and aesthetics is 
1/f noise and stochastic processes in nature generally. Given that the electroweak 
phase transition is stochastic, happening at random points in space governed by a 
temperature-dependent probability, it has a certain natural appeal for listeners. Follow- 
ing the discovery of 1/f noise in music and speech[4j, algorithms have been developed 
which can predict the "pleasantness" of musical and non-musical works[17|. Might 
their algorithm have some place in a program such as SuperCollider? 

Although SuperCollider is a powerful tool for sound synthesis, as a physics re- 
search project involving complex computations, Mathematica was the computer pro- 
gram of choice. Mathematica is common and well-established tool for computations 
but is clearly limited in its capacity for sound production. Likewise, there are many 
powerful and well-established computer programs in a variety of academic and non- 
academic sectors in which capacity for sound production, though likely beneficial, is 
limited or non-existent. The scope and speed of the evolution of sonification and au- 
ditory display would improve with the integration of sound production into these pro- 
grams. 
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